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Abstract 

We analyze the dispersion modes of hot plasma around Schwarzschild 
de-Sitter event horizon in the presence of Veselago medium. For this 
purpose, we apply ADM 3 + 1 formalism to develop GRMHD equa- 
tions for the Schwarzschild de-Sitter spacetime. Implementation of 
linear perturbation yields perturbed GRMHD equations that are used 
for the Fourier analysis of rotating (non-magnetized and magnetized) 
plasma. Wave analysis is described by the graphical representation 
of the wave vector, refractive index, change in refractive index, phase 
and group velocities. The outcome of this work supports the validity 
of Veselago medium. 
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1 Introduction 

Plasma is a distinct state of matter with a collection of free electric charge 
carriers which behave collectively and respond strongly to electromagnetic 
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fields. Debye shielding effect is the common example of collective effects of 
the plasma particles. The number density of positively charged particles re- 
mains almost the same as that of negatively charged particles to preserve the 
state of quasi-neutrality. Plasma is formed by the recombination of free elec- 
trons and ions, so it must have thermal energies to overcome the coupling of 
charged particles and can interact with electromagnetic fields. The nature of 
particle interactions distinguish between weakly and strongly ionized plasma. 
In weakly ionized plasma, charge-neutral interactions dominate the multiple 
Coulomb interactions, otherwise it is strongly ionized [1]. It is mentioned 
here that the most general plasma is hot and highly ionized [2]. 

Black holes are amongst the mysterious compact objects, depicting re- 
gion of space with strong gravitational field that nothing not even light can 
escape from it. According to the predictions of relativists, black holes are 
formed when matter is compressed to a point. These are commonly created 
from deaths of massive stars or by collapse of a supergiant star. Black holes 
contain the energy flux which generates a relatively large magnetic field [3]. 
The theory of magnetohydrodynamics (MHD) is developed for the study of 
plasma flow having both positive and negative charge carriers. This describes 
the study of plasma flow in the presence of magnetic field. The MHD theory 
along with the effects of gravity is said to be general relativistic magneto- 
hydrodynamics (GRMHD). The strong gravity of black hole perturbs the 
magnetospheric plasma. The theory of GRMHD is the most active field of 
study to examine the dynamics of magnetized plasma and impressions of 
black hole gravity. 

The de-Sitter black hole is a vacuum solution of the Einstein field equa- 
tions with a positive cosmological constant [4J. The Schwarzschild de-Sitter 
(SdS) spacetime expresses a black hole describing a patch of the de-Sitter 
spacetime. Inclusion of the positive cosmological constant in the Schwarzschild 
spacetime leads to SdS black hole which is non-rotating, thus plasma in mag- 
netosphere falls along the radial direction. Current observations indicate that 
our universe is accelerating, inclusion of positive cosmological constant leads 
to the expanding universe [5]-[Z]- Thus our universe approaches to de-Sitter 
universe in far future. 

Regge and Wheeler [8] and Zerilli [9] investigated the stability and grav- 
itational field of the Schwarzschild black hole with the help of non-spherical 
perturbations. Petterson [10] found that gravitational field is more strong 
near the surface of non-rotating black hole. Price [11] explored the dynamics 
of nearly spherical star by considering linear perturbation. Gleiser et al. [12] 
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discussed the stability of black holes by applying second order perturbations. 

The more appropriate approach towards General Relativity is to decom- 
pose metric field into layers of three-dimensional spacelike hypersurfaces and 
one-dimensional time, called as ADM 3 + 1 formalism, presented by Arnowitt, 
Deser and Misner (ADM) [13]. Thorne and Macdonald [H] examined that 
3 + 1 split is suitable approach for black hole theory. Smarr and York [TJ] 
used ADM formalism to explain spacetime kinematics numerically. Israel 
[16] studied event horizons in static vacuum and static electro- vacuum space- 
times. Zhang [T7] reformulated the laws of perfect GRMHD for the general 
line element and discussed the significant features of stationary symmetric 
GRMHD solutions. The same author [18] studied the dynamics of rotating 
black hole. Buzzi et al. [IS] investigated wave analysis of two fiuid plasma 
model around Schwarzschild event horizon. 

RezoUa et al. [20] considered the effects of cosmological constant and 
provided the dynamics of thick disks around SdS black hole. Myung [21] 
studied the entropy change for SdS black hole. Suneeta [22] presented the 
quasinormal modes for scalar field perturbations of SdS black hole. Ali and 
Rahman [23] investigated transverse wave propagation in two fiuid plasma for 
SdS black hole. Sharif and his collaborators [23]- [27] found dispersion modes 
of cold and isothermal plasmas for non-rotating as well as rotating black 
holes in the usual medium. They have also investigated the wave properties 
around Schwarzschild magnetosphere for the hot plasma. 

Artificial materials which exhibit unusual electromagnetic properties are 
termed as metamaterials. Veselago medium is the well-known metamate- 
rial named after Russian physicist Veselago [28j, which has negative elec- 
tric permittivity and magnetic permeability simultaneously. It has various 
names, double negative medium (DNM), backward wave medium (BWM), 
left-handed medium (LHM), negative refractive index medium (NIM) or neg- 
ative phase velocity medium (NPV). Many people [29]- [33] have explored the 
unusual properties of such a medium. Ziolkowski and Heyman [33] provided 
wave properties of NPV both analytically and numerically. Ramakrishna 
[55] discussed the importance of NIM for the perfect lensing. Sharif and 
Mukhtar [36] explored wave properties around Schwarzschild magnetosphere 
for isothermal as well as hot plasma in this unusual medium. 

In a recent paper (Sharif and Noureen, submitted), we have extended 
this study to SDS black hole for isothermal plasma. Here we investigate 
wave analysis of hot plasma for SdS black hole in NPV. For this purpose, 
the 3 + 1 GRMHD equations are considered to determine dispersion relations 
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with the help of Fourier analysis after insertion of linear perturbation to the 
GRMHD equations for both non-magnetized and magnetized backgrounds. 
The three dimensional plot of wave vector, refractive index and change in 
refractive index describes the wave properties of hot plasma. The paper is 
arranged as follows: In section 2, linearly perturbed 3 + 1 GRMHD equa- 
tions and their Fourier analysis is considered for hot plasma. The Fourier 
analyzed form of the GRMHD equations for rotating (non-magnetized and 
magnetized) plasmas are taken in sections 3 and 4 respectively to discuss the 
corresponding wave properties. Last section provides the conclusion. 

2 GRMHD Equations for Hot Plasma in a 
Veselogo Medium 

The general line element in ADM 3+1 formalism is represented as follows 

m 

ds'^ = -a^dt^ + r]ij{dx' + (3'dt){dx^ + (3^dt). (2.1) 

Here a denotes lapse function (ratio of FIDO proper time to universal time 
i.e., ^), a natural observer linked with this spacetime is called fiducial ob- 
server (FIDO), /3* describes three-dimensional shift vector and rjij = 
1,2,3) denote the components of three-dimensional spacelike hypersurfaces. 
The SdS spacetime in Rindler coordinates is given by 



ds'^ = -a^{z)dt^ + dx^ + dy"^ + dz^ , (2.2) 

where the directions z, y and x are analogous to the Schwarzschild coordi- 
nates r, (p and 9 respectively. The SdS black hole is non-rotating, thus the 
shift vector vanishes. From the comparison of Eqs. fl2.1l) and (12. 2p . we have 

a = a{z), (3 = 0, r],, = l{i=j). (2.3) 

Appendix provides the 3 + 1 GRMHD equations for the line element (12. 2p 
in a Veselago medium, i.e. Eqs. (]A5l) - (]A9p . The equation of state for hot 
plasma is [T7j 

, = (2.4) 

Po 

where the rest mass density, moving mass density, pressure and specific en- 
thalpy are denoted by po, p, p and p respectively. For hot plasma, the 
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specific enthalpy is not constant which shows that there is energy exchange 
between plasma and magnetic field of fiuid. The 3 + 1 GRMHD equations 
( flA5l) - flA9l) ) for hot plasma around the SdS black hole takes the form 

_ = -Vx (aVxB), (2.5) 

V.B = 0, (2.6) 

+ (p + P)y V.^ + (p + p)7V.(«V. V)V 

+ (p + p)V.(aV) = 0, (2.7) 

2 



+V.V) + 7V.(V.V)(p + p) - {^5., - l^B^B^j V;,V' 

= -{p + Ph'a, -p^^ + ^(yx B),V.(V X B) - ^(aB); 

^ [aS,),,B^--^[Bx{Vx(Vx(aVxB))}],, (2.^ 



4:7ra ' Aixa 

{-^ + V.V)(p + p)f - + 2(p + p)f (V.a) + {p + p) 
a at a ot 

-'^(V.V) - (V X B).(V X - - (V X B).(B x -) 

+ -^(V X B).(V X aB) = 0. (2.9) 

Ana 

We assume that for rotating background plasma fiow is in two dimensions, 
i.e., in xz-plane. Thus FIDO's measured velocity V and magnetic field B 
become 

Y = V{z)e^ + u{z)e,, B = S[A(^)e, + e,], (2.10) 

where B is an arbitrary constant. The following expression provides the 
relation between the quantities A, u and V [2] 



V= — + Am, (2.11) 

a 

where is constant of integration. The Lorentz factor, 7 = , , takes 
the form 

1 



Strong gravity of black hole perturbs the plasma flow. To determine the 
effect of black hole gravity on plasma flow, we apply linear perturbation. The 
flow variables (mass density p, pressure p, velocity V and magnetic field B) 
become 

p = + 5p = + pp, p = + 5p = + pp, 

V = V° + 5V = V° + V, B = B° + 5B = B° + 5b, (2.13) 

where p°, p, V° and B° represent the unperturbed quantities. The linearly 
perturbed quantities are denoted by Sp, Sp, (5V and 5B. We introduce the 
following dimensionless quantities p, p, Vx, Vz, and bz which correspond 
to the perturbed quantities 

p = p{t,z), p = p{t,z), V = SV = Vxit,z)ex + v^{t,z)ez, 

h = ^ = bx{t, z)ex + bzit, z)ez. (2.14) 

After the insertion of linear perturbations in the perfect GRMHD equations 
(Eqs. fl2.5l) - fl2.9p ) and using Eq. (12.141) . the component form of linearly per- 
turbed GRMHD equations become [36j 

1 dh 

— ^ - ubx,z = {ubx - Vbz -Vx + \Vz)V In a 
a at 

-{Vx,z - XVz,z - X'Vz + Vbz + Vbz,z - u'bx), (2.15) 

-§ = 0, (2.16) 
a at 

b,,, = 0, (2.17) 
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1 dp 1 dp . . o^^Advx . . , 2 

P-^+ P-jn + (P + Pn y{ — -KT + uv:,^,) + (p + p)7 u 
a ot a at a at 

1 uv 

+ (p + p)(l + 7V)i;,,, = -7Vp + p)[(1 + 27V2)\/' 

+272MyM>, + (p + - 272^2) (1 + 72^2) _ 

-27VF1/>;„ (2.18) 
+ |(p + p)7'(l + 7^2) + 1^1 + {{p + p)^^uV 

+^MpP + PP) {(1 + 7^'t^^)^' + I'^uVu'} + -f'^uV{p'p + pp' 
+P'P + pp') + [{P + P)l^u {(1 + 472^2)^^' + AW {I + 72^2^ I 

B'^ua' 

+ + 7'«(1 + 272^2) (p' + + + + 272^2) 

(1 + 27V')y' - 7W' + 27^(1 + 2^\^)uVu'} - l^(Aa)' 
+72^(1 + 2^\^){p' + p')]v, = 0, (2.19) 

(p + p)72(l + 'j^^) + —— }-^ + <{p + pHuV 



X 



An \ a dt \ An 
1 dv { X^B^ 1 

- + I (P + P)7' (1 + 7'«') + ^ j uv,,z +{{P + pHuV 

Ai?2 ^ XB"^ . r/ \ \/ /\ \/ / 

uVx,z H — -. — (1 + u )Ox^z + — {{^X) — a X + uX{ua 



An J ' An ' ' ' Ana 
+u'a)} bx + {pp + pp)j'^ {uz + uu'{l + 7^^) + -f'^u^VV'} 
+{1 + ^'^u^){p'p + pp') + -f^u^{p' p + pp') + [(p + p)7^ 

x{uW{l + 472^2) + 2V{az + uu'{l + 272^2))} _ 

Ana 

+27VV(p' + p')]vx + [(p + p)j^ {u'{l + + 472^2) 

+2^72(0, + (1 + 2^\'')VV')} + ^^{aXy + 272^1 

+7V)(p' + p')]^. = 0, (2.20) 
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1 dp 1 d'p 

Pttt + '^I'^uipp + pp)a2 + j'^u'ipp + pp) + 2(p + p)j^(uV' 

a ot 

+2uVa^ + u'V)vx + 2{p + p)-f'^{2'y'^uu' + a.-y^ + 2'fu^a^)v-, + 
2{p + p)^\Vv,,, + {p + p)j\l + 2^\^)v,^, - -—[{V^ + «')A^ 



dt ' ' dt ' ' (9t ^ 47ra 

X [\{\V + + (AF + u)v,^t - (A' + l)Vvx,t - (A' + 1)^^;.,*] 

+— (AA't;, - A'v, - + A'm6, - ^6,,, + u\K,z) = 0. (2.21) 

For the Fourier analysis of the hnearly perturbed GRMHD equations, we 
take the following harmonic spacetime dependence 



p{t, z) 


= cie~ 




-kz) 


p{t, z) 




Vz{t,z) 


= Cse" 


-i{uit~ 


-kz) 


Vx{t,z] 


1 = C4e~'(^*"^^) 


bzit,z) 


= Cse" 


-i(uit- 


-kz) 

1 


bx{t, z) 





(2.22) 

where u represents angular frequency and k denotes the z-component of the 
wave vector (0, 0, fc). The wave vector is helpful to determine refractive index 
and the properties of plasma about the event horizon. Frequency dependence 
effects on the wave propagating in a medium are related to dispersion. Using 
Eq.( l2.22| ) in Eqs.f l2.l5]) -f l2.2ip . we obtain their Fourier analyzed form 

C4(a' + ika) — C3 {{oiXy + ikaX } — c^{aV)' — CQ{{au)' — luj 
+ikua} = 0, (2.23) 

C5 = 0, 2.24 

a 

C5Lk = 0, (2.25) 

ci(^p) + C2(^p) + C3(p + p)[^-f^u + (1 + 7V)t/c 
a a a 

-(1 - 2j'u'){i + ^\^)- + 27Vyy'] + C4(p + ph% 



u a 



+Lku)V + u{l + 2YV')V' + 2Yu'Vu'] = 0, (2.26) 
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Ci[p7^m{(1 + 7^2)^' + + 7Vm(p' + ikp)] 

+C2b7^ii{(l + 7V^)y' + 7VW} + 7Vii(p' + iA;p)] 

+C3[(p + p)7'{(l + 27V)(1 + 27V)y' + + iku)jWu 

a 

-7W' + 27^(1 + 272^2)^^^'} + ^2y(l + 2-fhi^){p' + p') 
^\\ay + - ./cu)] + C4[(p + p)7M(l + 47^2) 



xuu' + 4VV'(1 + 7^2)} + (p + p)7^(l + 7V^)(^ + 

a 

+f «(1 + 27V^)(p' + ^>0 + ^ - ^(if^ _ ,ku)] 

Ana 47r a 

-C6- — [aW + a\l + -u^) + (1 + M2)/.A';a] = 0, (2.27) 

ci[p7^{a^ + (1 + 7^M^)W + j^u^VV'} + 7^M^(p' + tkp)] 
+C2W{a, + (1 + 7'w')W + 7WV^'} + (1 + 72^2) 

x{p' + ikp)] + C3[(p + p)7^{(l + 7V) (—^ + i/cw) 

a 

+u{l + 7^M^)(1 + 472^^) + 2u-f\a, + (1 + 272M2)\/y')} 

+2f u(l + 7V)(p' +pO + ^{Xay - - iku)] 

Ana 47r o; 

+C4[(p + Ph^'ii^ + iku)uV + uW{l + 47V) + 2V{a, 

a 

+(1 + 2j\^)uu')} + 2j%^V(p' + p') + ^(— - iku) 

Alt a 

4-770; 47rQ; 

+^(l + M2)tA;] = 0, (2.28) 



47r 

ci{(— ^7^ + ''^«7^ + 2uj^a, + 7^m')p + Vt^} + C2{( — (1 - 7^) 
a a 

+iku^'^ + 2^^ua, + -i^u')p + M7V} + C37^{(p' + p') + 2 

X (27W + + 2^^u^a^){p + p) + {l + 2-f^u^){p + p)ik + 

Ana 

x{\u - V)iu + aX'} + C4[2(p + p)-f^{{uV' + 2uVa, + u'V) + uVtk} 

B^ B^ 

H (V — uX)iUJ — aX'] + Cq[ {{ioou + iak)(Xu — V) — aX' 

Ana Ana 

+uViak}] = 0. (2.29) 



3 Rotating Non-Magnetized Background 



In the rotating non-magnetized plasma flow, we take B = = X and C5 = 
= C6 in the Fourier analyzed perturbed GRMHD equations f (l2:26|) - (l2:29|) ) 
[36]. 

3.1 Numerical Solutions 

To determine numerical solutions of the Fourier analyzed equations in rotat- 
ing non-magnetized background, we assume 



Specific enthalpy: = 
Time lapse: ot = 

Velocity components: u = V, x and z-components of velocity yield 
" = ^ = -^' 
Stiff fluid: p = p = f , 



vW2 



• Lorentz factor: 7 = , — , 

where is the SdS event horizon greater than that of the Schwarzschild 
event horizon and rn ^ 2M (l + ^ + ...j ^ C2M8km, l^C ^1.5 [23j. It 

is known that strong gravity of black hole plays an important role to perturb 
plasma present in its surroundings which causes the plasma wave propagation 
around black holes. The frequency dependence effects in wave propagation 
leads to the dispersion, i.e., the phenomenon in which phase or group velocity 
of the wave depends upon the frequency |38] . 

We consider the region 4 < 2; < 10 for wave analysis considering that 
event horizon is exactly at 2; = 0. The determinant of the coefficients of 
constants of the corresponding equations for the rotating non-magnetized 
plasma yields a complex dispersion relation. The real part of the determinant 
gives a quartic equation in k 

Ai{z)k^ + A2{z, u)k^ + Asiz, u)k^ + A^{z, u)k + A^iz, u) = (3.30) 

which yields two real and two imaginary roots. A cubic equation in k is 
obtained from the imaginary part of dispersion relation 

Bi{z)k^ + ^2(2, u)k'^ + Bs{z, u)k + B^{z, w) = (3.31) 
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which provides three real roots. The wave vector, refractive index, its change 
with respect to angular frequency, group velocity and phase velocity lead to 
the wave properties of the SdS black hole in the presence of Veselago medium. 
These are shown by graphical representation in Figures 1-5 expressing real 
values of k in Eqs.f l3.30|) and fl3.3ip . 

The dispersion is normal if phase velocity is greater than the group veloc- 
ity, otherwise anomalous [SS] or equivalently dispersion is normal if change in 
refractive index with respect to angular frequency is positive and anomalous 
otherwise. The information about energy exchange within the magnetosphere 
can be extracted only where the dispersion is normal. It is clear from figures 
that some waves move towards the event horizon and some move away from 
the horizon. The dispersion is normal in Figures 2, 3 and 5, while Figures 
1 and 4 indicate normal as well as anomalous at random points of the re- 
gion. The wave properties of non-magnetized hot plasma are illustrated in 
the tables I and II. 



Table I. Direction and refractive index of waves 



Fig. 


Direction of Waves 


Refractive Index (n) 


1 


Move towards the event horizon 


n < 1 and decreases in the region 
4<z<9,0<w<10 
with the decrease in z 


2 


Move away from the event horizon 


n < 1 and increases in the region 
4 < z < 9.8,0 < a; < 10 
with the decrease in z 


3 


Move towards the event horizon 


n < 1 and increases in the region 
A<z< 5.8,0 <uj <10 
with the decrease in z 


4 


Move towards the event horizon 


n < 1 and decreases in the region 
4:<z< 8.8,0 < w < 10 
with the decrease in z 


5 


Moves away from the event horizon 


n < 1 and increases in the region 
A<z< 5.8,0 <uj <10 
with the decrease in z 



The regions for anomalous and normal dispersion for figures exhibiting ran- 
dom dispersion are separated as follows: 

Table II. Regions of normal and anomalous dispersion 
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Fig. 


Normal dispersion 


Anomalous dispersion 


1 


4.6 < 2 < 10,0 < a; < 10 


A< z < 4.5,0.1 < w < 2.1 


4 


4.7 < z < 10,1 < w < 10 


4:< z < 4.5,0 < w < 3.5 



4 Plasma Flow With Rotating Magnetized 
Background 

Here plasma is considered to be rotating and magnetized. The flow is sup- 
posed to be in two dimensions, so the magnetic fleld and velocity of fluid lie in 
x2;-plane. Equations fl2.23p - fl2.29p represent the perturbed Fourier analyzed 
GRMHD equations for the rotating magnetized plasma. 

4.1 Numerical Solutions 

The assumptions for velocity, lapse function and speciflc enthalpy are the 
same as in the previous section. We put the following restrictions on magnetic 
field to find numerical solutions. 

• = 2 with u = V 

• putting = lm Eq.f l2TT]) so that A = 1 + 

Again our region of consideration is 4 < 2; < 10 and < w < 10. Equa- 
tions fi2.24p - fi2.25p yield C5 = 0. The dispersion relation, originated from the 
determinant of the Fourier analyzed form, gives the real part 

Ai{z)k^ + A2{z,u)k^ + A3{z,u)k'^ + A4{z,u)k + A5{z,u) = (4.1) 

which provides four roots, all are imaginary. The imaginary part of the 
dispersion relation 

Bi{z)k^ + B2{z, u)k'^ + B^{z, u)k^ + B^{z, u)k'^ + B^{z, uj)k 
+Be{z,uj) = (4.2) 

yields five roots of k, three real and two complex. The real roots indicate wave 
propagation in the region 4 < 2; < 10 and < a; < 10, shown in Figures 6-8. 
The waves are directed towards the event horizon given in Figure 6, move 
away from the event horizon as shown in Figures 7 and 8. The dispersion 
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Surface Plot of Have Vector Surface Plot of refractive Index 




Figure 1: Dispersion is normal and anomalous in the region. 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 2: Whole region exhibits the normal dispersion. 
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Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 3: Normal dispersion of waves is observed. 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 4: Random points of normal and anomalous dispersion are found in 
the region. 
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Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 5: Whole region admits normal dispersion. 



is normal as well as anomalous at random points in Figures 6 and 8, while 
Figure 7 indicates normal dispersion in the whole region. The tables III and 
IV summarize the results deduced from these figures. 

Table III. Direction and refractive index of waves 



Fig. 


Direction of Waves 


Refractive Index (n) 


6 


Move towards the event horizon 


n < 1 and decreases in the region 
6.95 <z< 7.3,0 < a; < 10 
with the decrease in z 


7 


Move away from the event horizon 


n < 1 and increases in the region 
4.3 < z < 4.7,0 < a; < 10 
with the decrease in z 


8 


Move away from the event horizon 


n < 1 and increases in the region 
4.4 < z < 4.9,0 < a; < 10 
with the decrease in z 



The regions of normal and anomalous dispersion for the roots exhibiting 
random dispersion can be separated as follows: 
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Surface Plot of wave Vector Surface Plot of Refractive Index 




Figure 6: Normal and anomalous dispersion of waves is observed. 

Surface Plot of wave Vector Surface Plot of Refractive Index 




Figure 7: Waves disperse normally in the whole region. 
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Surface Plot of wave Vector 



Surface Plot of 



Refractive Index 




Table IV. Regions of random dispersion 



Fig. 


Normal dispersion 


Anomalous dispersion 


6 


4<-z<7,2<a;<10 
7.9 < z < 10,3 < a; < 10 


1 <z< 7.8, 1 < a; < 4 


8 


h.l <z< 10,0.5 < a; < 10 


4 < z < 4.6,0.03 <uj< 0.35 
4.7 < 2 < 4.9,0.45 < w < 1.6 



5 Summary 

This paper is devoted to the study of dispersion modes for hot plasma sur- 
rounding SdS black hole in a Veselago medium. For the formulation of 
GRMHD equations in this unusual medium, we use the ADM 3 + 1 for- 
malism. The strong gravitational field of black hole disturbs the flow of 
surrounding plasma. To determine the impression of black hole gravity on 
plasma flow, we have applied linear perturbations to the GRMHD equations. 
Further, we have assumed that plasma flow is in two dimensions and evalu- 
ated the component form of linearly perturbed GRMHD equations. Finally, 
we have established dispersion relations for the rotating (non-magnetized and 
magnetized) plasma with the help of Fourier analysis techniques. 

In rotating non-magnetized plasma, waves move away from the event 
horizon shown in Figures 2 and 5 while waves are directed towards the event 
horizon given in Figures 1, 3 and 4. The dispersion is found to be normal 
in Figures 2, 3 and 5 in the whole region. It can be noticed that Figures 
1 and 4 exhibit normal as well as anomalous dispersion at random points. 
For rotating magnetized background, waves are directed towards the event 
horizon in Figure 6 while Figures 7 and 8 indicate that waves moves away 
from the event horizon. The dispersion is randomly distributed (normal as 
well as anomalous) shown in Figures 6 and 8 while Figure 7 admits normal 
dispersion in the whole region. 

In the usual medium, the refractive index is always greater than one but 
for Veselago medium its value must be less than one. The graphs of refractive 
index in all flgures show that region has refractive index less than one and 
increases or decreases in small regions. Also, the phase velocity is greater 
than group velocity for both non-magnetized and magnetized backgrounds. 
These are the prominent features of the Veselago medium which demonstrate 
the validity of this unusual medium for both rotating (non-magnetized and 
magnetized) background around SdS event horizon. 
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We have seen that for hot plasma surrounding the Schwarzschild black 
hole, most of the waves have anomalous and normal dispersion at random 
points while some waves exhibit anomalous dispersion in the whole region 
[36j . Moreover, most of the waves are directed away from the event horizon 
of the Schwarzschild black hole. However, here we have found that disper- 
sion is normal at most of the points. Thus it can be established that more 
information can be extracted from magnetosphere by the inclusion of pos- 
itive cosmological constant in non-rotating black hole. Also, the outcome 
of SdS black hole for isothermal plasma in a Veselago medium [37J shows 
wave propagation in different regions, i.e., —5 < z < —1, 0<a;<10 and 
l<z<5, 0<uj<10. For the isothermal plasma, some waves move 
towards event horizon and other move towards outer end of magnetosphere. 
This comparison shows that the wave properties can be better understood 
for the highly ionized, i.e., hot plasma around SdS magnetosphere. 

Appendix Al 

The Maxwell equations for such a medium are 




(Al) 
(A2) 



(A3) 



V X B 



(A4) 
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The GRMHD equations for the SdS spacetime in Rindler coordinates become 

_ = -V X («V X B), (A5) 
V.B = 0, (A6) 

-g + (aV.V)/9o + Po7'V.— + po7'V.(aV.V)V 

+PoV.(aV)=0, (A7) 



= -po/i7V -/>,i +;;^(V X B),V.(V x B) - -^{aB)\ 



S^' - -^[B X {V X (V X (aV X B))}]^, (A8) 
47rQ; 47rQ; 

+ V.V)( - - 1 + 2 W(V.a) + W(V.V) 
__L(v X B).(V X - —(V X B).(B X - V-) 

+-^(V X B).(V X aB) = 0. (A9) 
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